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Thermal  and  Oxidation  Response  of  UHTC  Leading  Edge  Samples 
Exposed  to  Simulated  Hypersonic  Flight  Conditions 
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Sharp  leading  edge  (LE)  samples  of  UHTC  (20  vol%SiC- 
HfB2)  and  SiC  were  exposed  to  simulated  hypersonic  flight 
conditions  using  a  direct-connect  scramjet  rig  and  their  ther¬ 
mal  and  oxidation  responses  measured.  The  measured  hack- 
wall  temperatures  and  scale  thicknesses  were  significantly 
smaller  than  might  he  expected  from  stagnation  temperatures 
at  the  LE.  Furthermore,  the  scale  that  formed  around  the  LE 
was  more  uniform  than  expected  from  the  steep  drop  in  cold 
wall  heat  flux  with  distance  from  the  tip.  These  results  were 
interpreted  and  rationalized  using  physics-hased  models.  An 
aerothermal  model  in  combination  with  an  oxidation  model 
accounted  for  the  observed  scale  thicknesses  at  the  tip  and 
their  slight  variation  with  distance.  The  scale  thicknesses  were 
similar  to  values  reported  for  exposures  in  furnaces  at  tem¬ 
peratures  calculated  for  the  tip,  but  less  than  those  reported 
in  arc  jet  tests.  The  formation  of  hafnon  (HfSi04)  and  the 
absence  of  external  glassy  layer  and  of  silica  in  the  outer 
portions  of  the  oxide  region  are  unique  to  scramjet  tested 
samples,  presumably  due  to  the  high  fluid  flow  (high  shear 
and  evaporation)  rates. 

I.  Introduction 

The  cowl  (LE)  and  fuel  injection  struts  of  a  scramjet 
engine  present  two  critical  materials  challenges  in  the 
development  of  reusable  hypersonic  flight  vehicles.  Both  are 
required  to  have,  and  retain  during  flight,  a  sharp  curvature 
(762  pm  or  30  mil  radius  in  current  design)  while  bearing  the 
extreme  conditions  of  hypersonic  flow.  The  conditions  experi¬ 
enced  by  a  sharp  body  under  hypersonic  flow  conditions 
have  been  the  subject  of  modeling  since  the  early  1950s.^“^  It 
is  widely  recognized  that  environmental  resistance  is  impor¬ 
tant  for  this  application. Key  parameters  that  represent 
these  conditions  with  respect  to  material  survivability  are 
heat  flux,  total  or  stagnation  temperature,  total  or  stagnation 
pressure,  dynamic  pressure,  fluid  velocity  near  the  material 
surface  behind  the  bow  shock,  fluid  composition,  the  degree 
of  dissociation  of  gaseous  elements  and  catalytic  recombina¬ 
tion  at  the  material’s  surface.  In  addition,  resistance  to  envi¬ 
ronmental  attack,  vibrations,  and  thermal  shock  will  also  be 
important.  However,  the  difficulty  in  reproducing  these  con¬ 
ditions  has  precluded  a  proper  evaluation  of  materials  for 
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this  application.  Most  of  the  experiments,  that  study  poten¬ 
tial  candidate  materials,  are  conducted  in  a  laboratory  fur¬ 
nace,  thus  assigning  the  total  temperature  as  the  key 
parameter  by  default. Alternately  it  is  possible  to  use  a 
laser  as  a  heat  source  where  candidate  materials  are  evalu¬ 
ated  by  comparing  performance  under  a  calculated  cold  wall 
heat  flux,  and  in  some  cases  with  imposed  fluid  flow.^^’^^  The 
arc-heater  or  arc-jet  test,  originally  developed  to  simulate 
atmospheric  reentry,  is  sometimes  used  to  evaluate  materials 
for  hypersonic  applications,  by  reproducing  a  calculated 
enthalpy. However,  the  heat  flux  in  these  tests  varies 
depending  on  the  material  surface-dependent  catalytic  coeffi¬ 
cient  for  recombination  of  gases.  Finally,  testing  using  real 
hypersonic  flights  is  very  expensive  and  thus  rare.^^ 

We  recently  reported  that  using  a  scramjet  rig  samples  can 
be  tested  under  conditions  that  simulate  to  a  reasonable 
accuracy  free-ffight  conditions. This  test  reproduces  the 
same  heat  flux,  total  temperature,  stagnation  pressure,  and 
fluid  velocity  behind  the  bow  shock  as  experienced  in  free 
flight  at  an  altitude  of  25  km.  Its  key  limitation  is  that  the 
gas  chemistry  is  different  from  air,  with  lower  oxygen  and 
higher  moisture  and  carbon-dioxide  levels.^  Among  the  can¬ 
didate  materials  that  are  expected  to  survive  these  conditions 
are  ultra-high  temperature  ceramics  (UHTC)  made  of  ZrB2 
or  HfB2  with  20  vol%  SiC,  which  have  been  shown  to  have 
good  oxidative  stability  up  to 

In  this  work,  experiments  were  conducted  that  explored 
the  thermal  and  oxidation  responses  of  LE  samples  made  of 
two  material  choices,  UHTC  and  SiC,  subjected  to  the  simu¬ 
lated  hypersonic  conditions  of  Mach  numbers  ranging  from 
4  to  7  at  an  altitude  of  25  km.  The  thermal  and  oxidation 
responses  of  the  samples  were  measured,  modeled,  and  com¬ 
pared  to  responses  in  furnace  tests  and  arc  jet  tests  reported 
in  the  literature  to  help  interpret/rationalize  all  of  the  obser¬ 
vations. 


II.  Experimental  Procedure 

The  LE  samples  were  wedge-shaped  with  a  762  pm  (30  mils) 
radius  of  curvature  and  a  12°  included  angle  [Fig.  1(a)],  as  per 
design  guidelines  for  current  scramjet  for  reusable  application. 
The  UHTC  sample  was  electro  discharge  machined  from  SPS 
(Spark  Plasma  Sintering)  processed  material  of  composition 
20  vol%  SiC-HfB2,  as  described  in  Ref.  15  The  SiC  sample 
was  sintered  SiC  (Hexaloy  SE,  St  Gobain,  Niagara  Falls,  NY) 
that  was  machined  (diamond  ground)  by  a  commercial  vendor 
(Bomas  Machine  Specialties,  Inc.,  Somerville,  MA).  The 
details  of  the  direct-connect  scramjet  rig  used  in  this  work  can 
be  found  in  Refs.  33,34  Briefly,  it  is  a  scramjet  engine  built  to 
study  supersonic  combustion  over  a  range  of  simulated  flight 
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conditions.  The  rig  is  supplied  with  compressed  air  that  is 
heated  using  a  combustion  heater  to  stagnation  temperatures 
simulating  flight  conditions  over  a  range  of  Mach  numbers 
from  3  to  7,  at  dynamic  pressures  up  to  95.8  kPa.  The  combus¬ 
tion  heater  is  fueled  with  natural  gas;  oxygen  is  added  in  suffi¬ 
cient  quantities  such  that  the  oxygen  mole  fraction  of  the  test 
gas  matches  that  of  clean  air.  A  supersonic  nozzle  then  acceler¬ 
ates  the  test  gas  to  the  proper  Mach  number  at  the  entrance  to 
the  scramjet  engine.  The  resulting  supersonic  flow  enters  the 
combustion  section  where  ethylene  is  used  as  fuel.  The  com¬ 
bustor  exhaust  gases  have  an  enthalpy  that  is  equivalent  to  free 
stream  hypersonic  flow.  It  must  be  clarified  that  the  gas  flow 
rate  itself  is  supersonic  and  not  hypersonic,  although  the 
enthalpy  content  is  the  same  as  hypersonic  flow  thus  simulat¬ 
ing  free  flight  at  Mach  numbers  up  to  7.  The  facility  includes  a 
variety  of  conventional  and  advanced  instrumentation  for 
accurate  documentation  of  combustor  inlet  and  exit  condi¬ 
tions.  The  rig  has  been  calibrated  to  obtain  profiles  of  pitot 
pressure,  total  temperature,  and  wall  static  pressure  distribu¬ 
tions  for  a  wide  range  of  inlet  conditions.  In  addition,  combus¬ 
tion  fluid  dynamics  modeling  of  the  rig  has  been  conducted  to 
obtain  predicted  conditions  within  the  rig,  and  found  to  be  in 
good  agreement  with  experiments. As  shown  in  our  prior 
work,^^  it  is  possible  to  suspend  leading  edge  samples  in  the 
combustor  exit  flow  path  of  the  scramjet  to  simulate  condi¬ 
tions  of  hypersonic  flight  at  an  altitude  of  25  km.  A  sample 
holder  assembly  that  holds  one  or  more  leading  edge  samples 
rigidly,  without  failing  after  multiple  exposures,  was  inserted 
into  the  combustor  exhaust  stream  through  an  access  port  in  a 
probe  housing  at  the  exit  of  the  combustor  section.  The  assem¬ 
bly  was  designed  and  built  using  an  iterative  procedure  that 
employed  modeling,  preliminary  trials,  and  experimental  trials 
as  described  in  Ref.  23 

The  design  used  in  this  work  consists  of  a  water-cooled 
inconel  tube  with  welded  sample  holder  subassemblies  all  of 


which  were  plasma  sprayed  with  a  zirconia  thermal  barrier 
coating  [Fig.  1(a)].  A  key  aspect  of  the  design  is  a  dovetail 
slot  in  the  subassemblies  to  accommodate  the  12°  included 
angle  of  the  leading  edge  (LE)  samples.  The  LE  sample  was 
instrumented  using  a  thermocouple  that  passed  through  a 
capillary  tube  from  the  back  of  the  water-cooled  inconel 
tube.  The  thermocouples  were  spot  welded  to  the  back  of  the 
inconel  tube  with  nichrome  strips  in  such  a  way  that  the  tip 
of  the  thermocouples  firmly  contacted  the  center  of  the  back 
of  the  LE  samples,  hereafter  termed  the  back  wall  of  the  LE 
sample.  No  attempt  was  made  to  attach  the  thermocouple  to 
the  sample  or  to  place  it  in  a  machined  recess.  Thus,  accu¬ 
racy  of  measurement  was  sacrificed  so  as  to  avoid  possible 
premature  failure  of  the  samples.  Pyrometric  measurement 
was  precluded  by  the  quartz  window  used  in  the  rig,  which 
has  high  absorption  in  the  frequency  range  of  interest. 

The  scramjet  rig  was  operated  at  different  conditions  by 
employing  facility  nozzles  with  different  Mach  numbers,  and 
by  varying  the  vitiator  stagnation  conditions  and  combus¬ 
tor’s  fuel/air  equivalence  ratio.  Parameters  of  importance  for 
aerothermal  effects  such  as  static  temperature,  static  pressure, 
and  total  temperature  were  collected  from  prior  calibrations 
of  the  rig.  Samples  of  UHTC  and  SiC  were  exposed  to  differ¬ 
ent  conditions.  The  total  temperature  was  used  to  calculate 
an  effective  Mach  number  for  free  flight  at  an  altitude  of 
25  km  using  the  equations  given  in  Ref.  1 

M=[//vWi'  (1) 

=  (2) 

where  U  is  the  fluid  velocity,  R  the  individual  gas  constant 
for  air  (286.9  J/kg-K),  Ti,  the  static  and  total/stagnation 


Fig.  1.  (a)  Schematic  sketch  of  the  Leading  Edge  sample,  and  the  sample  holder  used  in  this  work,  (b)  Schematic  of  thermal  loads  on  the 

leading  edge  sample  tested  in  the  scramjet.  (b)  The  variation  in  aerothermal  convective  heat  flux  with  location  along  the  surface  of  the  leading 
edge  shown  for  r^aii  =  0  (cold  wall)  for  effective  Mach  6.5  at  an  altitude  of  25  km. 
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temperatures,  respectively,  M  the  mach  number  and  y  the  ratio  of 
specific  heats  of  air;  the  static  temperature  as  a  function  of 
altitude  was  obtained  from  the  standard  atmospheric  table.^^ 
The  combustor  exit  gases  have  characteristics  that  reproduce 
nearly  the  same  heat  flux,  total  temperature,  fluid  velocity, 
total  pressure  behind  the  bow  shock  as  are  experienced 
behind  the  bow  shock  under  hypersonic  flow  at  an  altitude 
of  25  km,  as  shown  in  prior  work.^^  The  effective  Mach 
numbers  for  each  condition  of  tests  reported  in  this  work  are 
calculated  from  the  total  temperature  using  Eq.  (2).  The 
temperatures  of  the  back  wall  of  the  samples  were  continually 
monitored  during  the  runs  using  thermocouples.  The  samples 
were  removed  from  the  holder  and  examined  by  optical  macros- 
copy  (Nikon  DlOO,  Nikon  Inc.,  Melville,  NY),  scanning  electron 
microscopy  (FEI  Quanta  Mkll,  EEI  Inc.,  Hillsboro,  OR)  energy 
dispersive  spectroscopy,  EDS  (for  preliminary  analysis)  and  elec¬ 
tron  microprobe  analysis  (CAMECA  SXlOO,  Cameca  Instru¬ 
ments  Inc.,  Madison,  WI)  to  determine  the  thickness  and 
chemistry  of  the  scales  that  formed  on  them  as  a  function  of 
location  from  the  tip  of  the  sample.  X-ray  diffraction  (Rigaku 
D/MAX  B,  RU-200BH,  Rigaku  Americas,  Allison  Park,  PA) 
was  performed  on  the  scales  formed  on  the  slant  sides  of  the 
leading  edge  samples  to  identify  phases  present. 

III.  Models 


(1)  Thermal  Response 

Interpretation  of  the  experimentally  observed  data  required  a 
model  as  the  thermal  profile  of  the  samples  from  the  front  to 
back  could  not  be  measured  directly  during  the  run.  Due  to 
thermal  expansions  and  vibrations  present  during  the  run,  it 
was  not  possible  to  ensure  good  thermal  contact  between  the 
thermocouple  and  the  back  wall  of  the  samples.  Thus,  mod¬ 
eling  the  thermal  response  was  essential  to  interpret  the  data. 
The  model  was  also  helpful  in  evaluating  the  effectiveness  of 
the  rig  in  simulating  thermal  profiles  that  are  similar  to  what 
might  be  present  during  free  flight. 

The  thermal  responses  of  the  scramjet  rig  samples  as  well 
as  those  that  might  be  experienced  in  free  flight  were  com¬ 
puted  using  a  finite  element  model  (FEM,  Abaqus  6.11, 
Dassault  Systems,  Waltham,  MA)  that  incorporated  orienta¬ 
tion  and  distance-dependent  aerothermal  heating,  radiative 
heat  flux  from  the  sample  to  the  atmosphere,  and  conduction 
within  the  sample.  For  the  simulation  of  free-flight  condi¬ 
tions,  leading  edge  bodies  of  dimensions  50.8  mm  x  50.8  mm 
were  used  with  adiabatic  conditions  at  the  back  wall.  For  the 
scramjet  samples,  the  actual  dimensions  of  the  samples 
(12.7  mm  x  12.7  mm)  were  used,  along  with  appropriate 
thermal  boundary  conditions  for  the  back  end  of  the  sample, 
as  shown  in  Fig.  1(b).  As  specified  in  Fig.  1(b),  the  boundary 
conditions  at  the  gripped  end  assumed  that  the  thermal  con¬ 
ductivity  of  air  was  0.05  W/mK,  that  a  gap  of  15  micrometer 
exists  at  the  bolted  interface,  and  a  0.5  mm  gap  exists 
between  the  sample  and  the  dovetail  sample  holder,  based  on 
best  estimates  of  real  conditions. 

The  convective  aerothermal  heat  flux  at  the  leading  edge 
tip  was  computed  from  expressions  for  aerothermal  heat  flux 
given  in  the  work  of  Zoby.^^  Thus  the  aerothermal  heat  flux 
is  given  by 


W/m2  =  3.88  X  -  h^)  (3) 

(/!aw  -  /Jw),  kj/kg  =  Cp(rt  -  r*all)  (4) 


where  is  the  stagnation  heat  flux  on  the  material,  pa  is 
the  total  pressure  behind  the  bow  shock,  r  the  radius  of  cur¬ 
vature  of  the  cylindrical  sample,  /^aw  the  adiabatic  wall 
enthalpy,  the  enthalpy  of  air  calculated  at  the  material 
wall  temperature,  Cp  the  specific  heat  of  air,  f  the  total  tem¬ 


perature,  and  Twaii  the  material  wall  temperature.  For  Mach 
numbers  of  up  to  7,  the  degree  of  dissociation  is  negligible 
(<1%)  for  both  O2  and  N2,  as  calculated  using  the  expres¬ 
sions  given  by  Scala  and  Gilbert^^  as  shown  in  Ref.  23  The 
stagnation  or  total  pressure  behind  the  bow  shock,  p^,  is 
given  by  the  Raleigh  pitot  tube  formula: 


Pt2  =  P\ 


f  (y+1)M^  \ 
\2yM^  -  (Y-  1)/ 


(5) 


where  pi  is  the  free  stream  (static)  pressure  obtained  from 
prior  calibrations  of  the  rig,  as  detailed  in  Ref.  23  Note  that 
the  heat  flux  varies  with  wall  temperature  and  this  is 
accounted  for  in  the  FEM  model.  The  heat  flux  also  varies 
with  distance  away  from  the  tip  due  to  a  change  in  the  orien¬ 
tation  of  the  surface  normal  with  flow,  and  is  given  by: 


G(0)  =  G,Cos(0)  (6) 

where  0  is  the  orientation  of  the  plane  normal  relative  to  the 
hypersonic  flow  direction,  which  varies  from  0°  at  the  tip  of 
the  EE  and  reaching  0c  (=84°)  as  the  curvature  meets  the 
slant  surface.  With  further  distance  along  the  wedge  surface 
at  constant,  0c,  the  heat  flux  drops  as 

Q{x)  =  Q^Cos(Q,)/V^  (7) 

where  v  is  the  distance  from  the  point  where  the  EE  curva¬ 
ture  meets  the  slant  surface  (the  tangent  point).  A  plot  of  this 
variation  in  heat  flux  is  shown  in  Fig.  1(c)  for  =  0  (cold 

wall)  for  an  effective  free  stream  Mach  number  of  6.5  at  an 

altitude  of  25  km.  It  must  be  noted  that  while  this  distribu¬ 
tion  in  flux  includes  the  correct  dependencies,  it  is  an  approx¬ 
imation  of  the  actual  distribution,  especially  the  transition 
from  angular  dependence  [Eq.  (6)]  to  distance  dependence 
[Eq.  (7)].  The  radiative  heat  flux  from  the  surface  to  ambient 
is  given  by  the  Stefan-Boltzmann  law  as  shown  in  Fig.  1(b). 
It  varies  with  material  surface  temperature  and  ambient  sta¬ 
tic  temperature.  The  transient  conduction  within  the  sample 
is  solved  using  the  finite  element  model.  The  thermal  conduc¬ 
tivity  of  the  interface  between  the  sample  and  sample  holder 
was  estimated  using  the  conductivity  of  air  (0.05  W/mK)  and 
an  estimated  gap  thickness,  as  shown  in  Fig.  1(b).  The  back 
wall  heat  flux  was  estimated  using  the  conductivity  of  air 
across  a  static  air  gap  and  radiation  to  the  water-cooled  inc¬ 
onel  maintained  at  an  estimated  temperature  of  450  K.  The 
results  at  the  tip  and  hotter  regions  of  the  sample  were  not 
sensitive  to  the  back  wall  heat  conduction  estimates.  For  the 
simulation  of  free-flight  conditions,  the  convective  aerother¬ 
mal  heat  flux  was  applied  to  the  entire  slant  surface,  and  an 
adiabatic  condition  was  assumed  for  the  back  wall. 

The  conditions  in  the  rig  such  as  total  temperature,  static 
temperature,  fluid  flow  velocity,  static  pressure,  and  the  spe¬ 
cific  heat  ratio  of  the  gas,  were  obtained  from  prior  calibra¬ 
tions  of  the  scramjet  rig.  Note  that  the  radiation  to  the 
ambient  in  the  scramjet  rig  is  dependent  on  the  static  temper¬ 
ature  present  in  the  rig,  but  in  the  free-flight  simulation  it 
depends  on  the  ambient  temperature  at  an  altitude  of  25  km, 
taken  from  standard  atmospheric  tables. 


(2)  Oxidation  Response 

The  oxidation  response  was  evaluated  using  electron  micros¬ 
copy,  energy  dispersive  spectroscopy  (EDS),  electron  probe 
microanalysis  (EPMA),  and  X-ray  diffraction  (XRD),  as 
mentioned  earlier.  The  resulting  data  were  interpreted  by 
comparing  the  scale  thicknesses  to  those  predicted  by  an  oxi¬ 
dation  model  developed  based  on  oxidation  experiments  con¬ 
ducted  in  furnace  tests.  The  model  used  in  this  work  is 
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detailed  elsewhere. The  model  was  found  to  predict  results 
from  furnace  tests  reasonably  well  at  lower  temperatures  (up 
to  1773  K).  It  was  found  to  underpredict  slightly  the  oxida¬ 
tion  rates  at  higher  temperatures,  partly  due  to  possible  con¬ 
tamination  from  the  high-temperature  furnaces  used. 
Furthermore,  the  model  was  found  to  significantly  underesti¬ 
mate  the  oxidation  rates  in  arc  jet  tests,  where  the  enthalpy 
is  carried  by  dissociated  gases  which  recombine  at  the  mate¬ 
rial  surface.  The  model  accounts  for  the  high  fluid  velocity  of 
the  ambient  present  behind  the  bow  shock  through  boundary 
layer  evaporation.  The  model  does  not  include  the  shear 
stresses  that  might  strip  the  liquid  layer,  but  under  the  condi¬ 
tions  of  the  test  the  model  predicts  the  evaporative  flux  to  be 
so  high  that  an  external  glassy  layer  barely  forms.  It  is  thus 
applicable  for  the  scramjet  conditions  used  in  this  work.  As 
mentioned  earlier,  gas  dissociation  and  recombination  is  not 
significant  at  Mach  numbers  below  7,  based  on  calculations 
using  equations  given  by  Scala  and  Gilbert. The  model  pre¬ 
dictions  for  oxide  layer  thickness  were  compared  with  those 
obtained  using  the  scramjet  rig  tests.  The  predictions  were 
based  on  thermal  responses  calculated  for  the  tip  of  the  lead¬ 
ing  edge,  using  the  FEM  model  detailed  earlier. 

IV.  Results 

(1)  Experimental  Results 

A  summary  of  all  the  runs  where  LE  samples  survived  are 
shown  in  Table  I.  The  conditions  of  the  run  were  obtained 
from  the  data  collected  during  the  run  and  using  prior  cali¬ 
brations  of  the  rig.  These  data  were  converted  using  appro¬ 
priate  equations  to  calculate  the  cold  wall  heat  flux  on  the 
leading  edge  of  the  sample  which  had  a  762  micrometer 
radius  of  curvature. (see  Ref.  5  for  details)  The  dynamic  pres¬ 
sure,  total  pressure,  and  fluid  velocity  behind  the  shock  were 
other  parameters  used  for  this  evaluation.  By  comparing  to 
the  total  temperatures  during  free  flight  at  an  altitude  of 
25  km,  the  conditions  in  the  rig  were  converted  to  an  effec¬ 
tive  Mach  number.  Figure  2  shows  a  plot  of  sample  survival 
conditions  (total  temperature  and  effective  Mach  numbers) 
for  the  UHTC  and  SiC  samples  tested.  UHTC  samples  were 
found  to  survive  more  often  than  SiC  samples.  The  inset 
shows  video  captures  taken  before,  during,  and  after  a  typi¬ 
cal  run.  For  comparison,  it  must  be  noted  that  a  graphite 
sample  tested  during  one  of  the  runs  (run  ID#  10-6  in 
Table  I)  burned  at  a  Mach  number  of  4.7  (using  vitiator  but 
no  combustion)  after  nearly  5  min  of  hold. 

Figures  3(a)  and  (b)  show  SEM  images  of  the  oxide  scales 
formed  on  the  UHTC  samples  for  a  typical  run  (9.3  min  at  a 
total  temperature  of  1550  K,  followed  by  13.3  min  at  total 
temperature  of  1850  K,  runs  #10-1 +#10-2)  as  a  function  of 
the  location  of  the  leading  edge  surface  with  respect  to  the 
tip.  Figure  3(c)  shows  the  X-ray  diffraction  pattern  obtained 
on  one  of  the  slant  sides  of  the  wedged  sample.  It  shows  the 
presence  of  both  hafnia  (Hf02)  and  hafnon  (HfSi04)  in  the 
sample.  Figure  4  shows  a  higher  magnification  of  the  scale 
formed  on  a  leading  edge  sample  after  a  different  run 
(4.33  min  at  a  total  temperature  of  1900  K  followed  by 
1.33  min  at  2370  K).  The  figure  includes  EPMA  results  col¬ 
lected  on  two  different  regions  of  the  scale  showing  regions 
of  Hf-Si-O  and  Hf-O.  Different  filters  were  used  to  identify 
different  elements  as  Si  and  Hf  peaks  overlap  on  EDS  maps. 
From  microprobe  (EPMA)  analyses  of  the  different  regions, 
combined  with  the  XRD  data,  the  scale  was  found  to  be 
comprised  mostly  of  hafnia,  with  an  outer  layer  containing 
hafnon.  Thus,  from  the  SEM  images,  EPMA  of  different 
regions,  and  XRD  of  the  scales,  the  oxide  scales  were  found 
to  contain  hafnia  and  hafnon,  with  limited  presence  of  silica, 
but  no  external  glassy  layer.  There  was  no  SiC-depleted 
region  beneath  the  oxide  scale.  In  contrast  to  UHTC  sam¬ 
ples,  no  scale  was  found  on  the  SiC  samples  that  survived. 
This  could  result  from  either  high  shear  stresses  or  high 
evaporation  rates  both  arising  from  the  high  flow  rate  condi¬ 


tions.  This  is  consistent  with  the  observation  that  the  sur¬ 
faces  of  the  scales  on  UHTC  samples  were  free  of  silica  and 
the  inside  of  the  scale  was  lean  in  silica.  The  formation  of 
hafnon,  the  lack  of  an  external  glassy  layer  and  the  lack  of 
silica  in  the  outer  regions  of  the  oxide  scale  are  features 
unique  to  scramjet  testing,  presumably  due  to  the  high  fluid 
flow  conditions,  which  can  cause  enhanced  loss  of  silica  from 
the  surface. 


(2)  Model  Results 

Figures  5(a)  and  (b)  show  the  material  thermal  properties 
used  in  the  model.  The  specific  heat  data  for  both  SiC  and 
UHTC  are  well-established. The  thermal  conductivity  for 
sintered  SiC  is  also  known  with  reasonable  certainty. The 
emissivity  of  UHTC  was  taken  to  be  0.75  based  on  the  work 
by  Meng  et  al?^  and  that  of  SiC  to  be  0.9.  The  model  results 
had  a  very  weak  dependence  on  emissivity,  but  a  strong 
dependence  on  specific  heat  and  thermal  conductivity.  The 
specific  heat  is  well-known,  but  a  large  variation  in  the  ther¬ 
mal  conductivity  for  UHTC  has  been  noted  in  recent  work- 
s.(see  e.g..  Ref.  18)  Hence  the  thermal  diffusivity  of  the 
sample  used  in  this  study  was  measured  using  laser  flash  up 
to  a  temperature  of  1200°C.  As  shown  in  Fig.  5(b),  this 
datum  converted  to  thermal  conductivity  (from  specific  heat 
and  density)  was  found  to  lie  very  close  to  an  average  curve 
fit  to  available  literature  data  in  Ref.  18  Figure  5(c)  shows 
the  typical  FEM  model  results  for  the  calculated  thermal 
response  of  the  UHTC  LE  samples  within  the  scramjet  rig 
and  that  of  a  UHTC  LE  in  free  flight.  The  scramjet  sample 
was  simulated  using  the  actual  dimensions  of  the  sample 
(30  mil  curvature  at  the  tip,  12°  included  angle, 
12.7  mm  x  12.7  mm)  and  sample  holder  assembly.  The  ther¬ 
mal  response  of  an  UHTC  LE  under  free-flight  conditions 
was  calculated  using  the  same  tip  curvature  and  included 
angle  of  the  wedge,  but  with  a  significantly  larger  sample  size 
of  50.8  mm  X  50.8  mm.  The  back  wall  of  the  sample  was 
taken  to  be  adiabatic.  From  the  plot  of  temperature  versus 
distance  from  the  leading  edge  tip  shown  in  Fig.  5(d),  it  is 
seen  that  the  sample  in  the  scramjet  rig  has  almost  the  same 
profile  as  that  in  free  flight.  The  LE  sample  is  cooler  only  in 
the  region  where  the  sample  holder  grips  the  sample,  as  is  to 
be  expected.  For  a  free  flight  of  Mach  7,  the  leading  edge 
reaches  a  temperature  of  -1850  K,  and  decreases  with  dis¬ 
tance  at  a  rate  of  -40  K/mm. 

Figure  6  shows  the  various  recorded  and  predicted  tem¬ 
peratures  in  the  rig  as  a  function  of  effective  Mach  number 
for  the  UHTC  and  SiC  samples.  The  plot  compares  the  mea¬ 
sured  back  wall  temperatures  measured  during  the  run,  the 
predicted  LE  tip  temperature  from  the  model,  and  the  pre¬ 
dicted  back  wall  temperature  from  the  model.  The  measured 
temperatures  are  slightly  below  the  predicted  temperatures, 
presumably  due  to  the  lack  of  proper  contact  and  effects 
from  thermal  expansion  and  vibrations,  as  explained  earlier. 
Note  that  the  LE  tip  temperatures  for  SiC  are  predicted  to 
be  slightly  higher  than  those  for  UHTC,  and  the  back  wall 
temperatures  lower,  quite  likely  from  the  lower  thermal  con¬ 
ductivity  of  SiC  at  high  temperatures. 

Figure  7  shows  a  comparison  of  the  oxide  thicknesses 
measured  using  scanning  electron  microscopy  on  the  LE 
samples  for  different  runs  of  UHTC,  with  the  curves  pre¬ 
dicted  by  the  model.  The  model  predictions  included  the 
effects  of  increased  total  pressures  behind  the  shock  wave, 
and  the  high  fluid  velocity  behind  the  shock  wave  that  the 
material  will  experience.  Data  from  furnace  or  arc  jets  were 
not  available  for  the  exact  temperatures  experienced  by  the 
leading  edge,  but  data  available  for  temperatures  that  are 
close  are  shown  for  comparison  purposes.  In  Fig.  7(b),  oxide 
thicknesses  measured  on  samples  tested  in  a  furnace  near  the 
same  temperature  is  included.  Figure  7(d)  shows  a  compari¬ 
son  of  UHTC  samples  tested  in  a  furnace  and  those  tested  in 
an  arc  jet  at  the  same  temperature.  The  model  predicts  that 
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Table  I.  A  List  of  Runs  with  Parameters  and  Calculated  Conditions,  Under  Which  Leading  Edge  Samples  Survived  Exposures  in  the  Scramjet  Rig 
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-  UHTCl*  sample  used  in  run#  11-4  was  reused  in  run#  11-5. 


912 


Journal  of  the  American  Ceramic  Society — Parthasarathy  et  al. 


VoL  96,  No.  3 


the  evaporation  rate  of  external  glass  will  be  sufficiently  fast 
to  prevent  a  stable  external  glassy  film,  as  experimentally 
observed. 

V.  Discussion 

This  work  presents  the  thermal  and  oxidation  responses  of 
UHTC  and  SiC  leading  edge  samples  exposed  to  hypersonic  flow 
conditions  which  simulated  effective  free  flight  Mach  numbers  of 
6-7  at  an  altitude  of  25  km.  The  temperatures  reached  during 
the  run  at  the  back  wall  of  the  sample,  which  is  12.7  mm  (0.5") 
from  the  leading  edge,  were  measured  using  thermocouples 
which  firmly  touched  the  back  wall  of  the  samples  at  room  tem- 


5.25  5.75  6.25  6.75  7.25 


Effective  Mach  Number 


Fig.  2.  A  plot  of  conditions  under  which  UHTC  and  SiC  leading 
edge  samples  were  exposed.  The  green  and  red  symbols  denote 
conditions  where  the  samples  survived  or  were  lost  respectively,  with 
circles  representing  UHTC  and  squares  representing  SiC.  UHTC 
samples  were  found  to  have  a  better  survival  rate  than  SiC  samples. 
The  inset  shows  video  captures  before,  during  and  after  a  typical  run. 
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perature.  The  scale  thicknesses  formed  on  the  samples  were  char¬ 
acterized  using  SEM,  EDS,  EPMA,  and  XRD.  The  results  were 
interpreted  using  a  combination  of  aerothermal,  thermal,  and 
oxidation  modeling.  The  following  is  a  discussion  of  the  results. 

The  scramjet  rig  conditions  vary  from  run  to  run,  but  all 
the  parameters  of  importance  for  aerothermal  calculations  are 
available  either  through  live  recordings  of  data  or  prior  cali¬ 
brations  to  control  parameters.  The  total  temperature  was 
used  to  calculate  an  effective  Mach  number  for  free  flight  at 
this  altitude.  Erom  Eig.  2,  it  is  clear  that  we  can  obtain  a 
range  of  Mach  numbers  spanning  from  5  to  7  under  combus¬ 
tion  conditions  which  the  UHTC  samples  are  capable  of  sur¬ 
viving.  The  thermal  modeling  results  shown  in  Eig.  5  illustrate 
that  the  thermal  profile  experienced  by  the  12.7  mm  (0.5") 
leading  edge  sample  in  the  rig  are  very  close  to  that  for  a 
much  larger  50.8  mm  (2")  sample  under  free-flight  conditions. 

When  combined,  the  experiments  and  modeling  show  that 
the  tip  of  the  UHTC  LE  sample  reaches  temperatures  far 
below  that  of  the  total  temperature.  As  shown  in  Eig.  6,  the 
measured  and  predicted  back  wall  temperatures  are  much 
lower  than  the  total  temperature  and  the  LE  tip  temperatures 
are  predicted  to  be  in  between  the  two.  Eor  example,  it  is 
seen  that  for  UHTC,  the  LE  tip  reaches  a  temperature  of 
1850  K,  when  the  total  temperature  is  2380  K  and  the  static 
temperature  in  the  rig  is  1450  K.  The  same  profile  is  reached 
in  a  larger  LE  of  UHTC  under  free  flight  conditions  of  Mach 
7  at  an  altitude  of  25  km  with  a  static  ambient  temperature 
of  220  K  as  seen  from  Eig.  5.  Thus,  for  reusable  hypersonic 
vehicles  the  maximum  temperature  reached  by  a  thermal  pro¬ 
tection  system  using  UHTCs  will  be  1850  K  for  up  to  Mach 
7.  Two  parameters  of  interest,  temperature  and  heat  flux,  are 
often  used  in  test  methods  to  evaluate  candidate  materials. 
Eigure  8  gives  suggested  values  to  be  used  in  such  tests  for 
various  free-flight  conditions.  The  total  temperature  and  tem¬ 
perature  reached  at  the  tip  of  the  LE  are  shown  in  Eig.  8(a). 
The  initial  cold  wall  heat  flux  and  steady-state  heat  flux  are 


(b) 


2.0  mm  from  tip  0.9  mm  from  tip 


Fig.  3.  (a)  Oxidation  scale  at  the  tip  of  a  UHTC  leading  edge  sample  after  exposure  in  the  scramjet  rig  for  9.3  min  at  a  total  temperature  of 
1550  K,  followed  by  13.3  min  at  total  temperature  of  1850  K.  (b)  oxidation  scale  formed  over  the  slant  edges  of  the  12°  wedge  away  from  the 
leading  edge  tip.  There  was  no  evidence  of  external  glassy  layer  or  internal  depletion  layer,  which  have  been  reported  to  form  in  furnace 
experiments,  (c)  X-ray  diffraction  pattern  obtained  from  the  slant  sides  showing  presence  of  hafnia  and  hafnon  in  the  scale. 
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Fig.  4.  Oxidation  scale  at  the  tip  of  an  UHTC  leading  edge  sample  after  exposure  in  the  scramjet  rig  for  4.33  min  at  a  total  temperature  of 
1900  K  followed  by  1.33  min  at  2370  K.  There  was  no  evidence  of  an  external  glassy  layer  or  internal  depletion  layer.  EPMA  analysis  of  the 
scale  showed  the  presence  of  both  hafnia  and  hafnon  in  the  oxide  scale.  Note  that  the  scale  used  for  Si  in  the  EPMA  corresponding  to  hafnia 
region  is  much  smaller  than  that  shown  for  the  hafnon  region.  X-ray  diffraction  of  the  scale  also  showed  the  presence  of  both  hafnia  and  hafnon 
in  the  scale. 


(a) 


(b) 


(c) 


(d) 


Fig.  5.  (a,b)  The  temperature-dependent  thermal  properties  used  in 

the  model  were  obtained  from  literature  and  experimental  data  on 
UHTC  used  in  this  work.  (c,d)  The  calculated  thermal  response  of 
the  leading  edge  in  the  scramjet  rig  is  shown  compared  with  the 
calculated  response  of  a  larger  sample  in  free  flight  using  adiabatic 
conditions  for  the  back  wall.  It  is  seen  that  the  thermal  response  of 
the  sample  in  the  scramjet  rig  follows  closely  that  in  free  flight 
except  for  the  portion  of  the  sample  that  is  held  in  the  cooled 
dovetail  region. 


shown  in  Fig.  8(b).  The  variation  in  material  temperature 
with  distance  from  the  tip  is  shown  in  Fig.  8(c).  The  scale 
thickness  formed  over  the  LE  is  shown  in  Fig.  8(d).  Note 


Effective  Mach# 

Fig.  6.  The  various  recorded  and  predicted  temperatures  for  all  the 
runs  are  shown  plotted  as  a  function  of  effective  Mach  number  of 
the  simulation  conditions  used  in  the  scramjet  rig  for  UHTC  samples 
and  SiC  samples.  The  higher  thermal  conductivity  of  UHTC  results 
in  lower  EE  tip  temperature  than  the  SiC  sanple. 

that  the  scale  thickness  does  not  vary  significantly  with 
distance  from  the  tip  for  short  exposures  (10  min)  consistent 
with  the  small  variation  in  scale  thicknesses  observed  on  the 
LE  samples  as  shown  in  Figs.  3  and  4. 

The  comparisons  of  the  oxidation  model  predictions  with 
the  experimental  data  from  samples  tested  in  the  scramjet 
show  reasonable  agreement,  but  the  experimental  values  are 
consistently  higher  than  the  predicted  values  by  a  small  mar¬ 
gin.  This  is  true  for  furnace  tested  samples  as  well,  but  only 
at  higher  temperatures  where  contamination  from  the  fur¬ 
nace  was  thought  responsible  for  the  discrepancy.^^  The  arc 
jet  tests  typically  show  a  larger  scale  thickness  than  furnace 
tests  as  shown  by  other  investigators  and  as  seen  from 
Fig. 8(d)  and  pointed  out  in  more  detail  in  Ref.  32  One  possi¬ 
ble  reason  for  the  observed  higher  oxidation  rates  in  the 
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Fig.  7.  The  oxidation  scale  thicknesses  measured  on  the  LE  samples  after  exposures  in  the  scramjet  rig  are  shown  plotted  with  model 
predictions,  furnace  oxidation  data,  and  arc  jet  oxidation  data.  The  scales  on  the  LE  samples  are  slightly  thicker  than  model  predictions  using 
the  predicted  material  temperatures,  but  not  very  different  from  the  trends  in  how  the  model  compares  with  the  furnace  data  at  high 
temperatures.  On  the  contrary,  the  arcjet  data  is  considerably  larger  than  both. 


Mach  # 


Mach  # 


Distance,  m  Distance,  m 

Fig.  8.  Calculated  thermal  response  of  UHTC  as  leading  edge,  plotted  as  a  function  of  effective  free  stream  Mach  number  at  an  altitude  of  25  km. 
(a)  Total  and  material  tip  temperatures,  (b)  Cold  wall  and  steady  state  heat  flux,  (c)  temperature  profiles  with  distance  from  the  tip  of  the  LE,  and 
(d)  predicted  thickness  of  oxide  scale  as  a  function  of  distance  from  the  LE  tip  for  2  min,  10  min  and  1  h  at  a  free-flight  simulation  of  Mach  7.  Note 
that  at  short  durations,  the  scale  thickness  does  not  vary  significantly  with  distance,  rationalizing  the  experimental  observations  (Fig.  3  and  4). 


scramjet  compared  to  the  model  prediction  even  at  lower 
temperatures  is  that  the  chemistry  of  the  scramjet  fluid 
includes  significant  moisture.  The  other  possible  reason  could 


arise  from  the  observation  of  the  formation  of  hafnon  in  the 
scale  which  would  result  in  volume  shrinkage  making  the 
scale  open  up  the  pore  fraction.  Finally  much  less  glassy 
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phase  was  found  on  the  scale  compared  to  the  samples  tested 
in  the  furnace.  It  is  believed  that  the  hypersonic  flow  is  able 
to  evaporate  or  pull  the  glassy  material  out  of  the  scale  thus 
resulting  in  a  porous  oxide  scale  with  no  glassy  phase  filling 
it,  at  least  at  the  outer  most  regions  of  the  scale.  Further 
characterization  of  the  scales  will  be  required  to  understand 
the  mechanisms  contributing  to  the  slightly  enhanced  oxida¬ 
tion  under  the  scramjet  and  arc  jet  conditions.  Despite  these 
limitations,  the  model  captures  the  dependencies  well  and 
predicts  oxidation  rates  within  reasonable  error  for  use  in 
engineering  design. 

VI.  Summary 

Sharp  leading  edge  samples  of  UHTC  (20  vol%SiC-HfB2) 
and  SiC  were  exposed  to  simulated  hypersonic  free-fiight 
conditions  using  the  post  combustion  section  of  a  direct-con¬ 
nect  scramjet  rig,  which  served  effectively  as  a  hypersonic 
wind  tunnel.  The  UHTC  samples  survived  more  often  than 
SiC  samples  and  withstood  exposures  for  up  to  15  min  under 
conditions  that  simulated  free  flight  at  an  altitude  of  25  km 
at  Mach  6.25,  and  survived  a  80  s  exposure  at  Mach  7.  The 
measured  back  wall  temperatures  were  far  below  the  total 
temperatures  and  the  oxidation  scales  were  thinner  than  what 
might  be  expected  to  form  at  total  temperatures.  These 
results  were  rationalized  using  a  combination  of  thermal  and 
oxidation  models.  The  tip  of  the  LE  samples  reached  temper¬ 
atures  calculated  to  be  350-500  K  below  the  total  tempera¬ 
tures  and  maintained  a  steep  gradient  of  40  K/mm  away 
from  the  tip,  at  steady  state.  Similarly  the  steady-state  heat 
fluxes  were  15%-20%  below  the  cold  wall  heat  fluxes.  The 
thermal  model  predicted  that  for  the  same  conditions,  the  tip 
temperature  for  UHTC  was  lower  than  that  for  SiC  due  to 
the  higher  thermal  conductivity  of  UHTC  at  temperature. 
The  scale  thicknesses  on  UHTCs  were  within  reasonable  lim¬ 
its  of  predictions  from  the  oxidation  model.  Possible  reasons 
for  slightly  enhanced  oxidation  in  the  scramjet  compared  to 
the  model  are  likely  due  to  the  formation  of  hafnon  in  the 
scale  and  the  impurities  in  the  scramjet  environment  includ¬ 
ing  higher  moisture  content.  The  scale  thicknesses  were  simi¬ 
lar  to  values  reported  for  exposures  in  furnaces  but  less  than 
those  reported  in  arc  jet  tests.  The  absence  of  an  external 
glassy  layer  and  absence  of  silica  in  the  outer  portions  of  the 
oxide  region  are  unique  to  samples  tested  in  the  scramjet, 
presumably  due  to  the  high  fluid  flow  rates. 
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